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CHAPTER  1  INTRODUCTION 

1.1  INTRODUCTION 

In  the  Air  Force  weapon-delivery  and  defense,  advanced  MHD  electric  power 
generators  play  an  Important  role.  In  supersonic  generators  using  high  magnetic 
field  strength,  designed  for  high  power-density*  the  Influence  of  the  wall 
roughness  Is  felt  rather  acutely  In  the  core  flow  because  of  propagation  of 
disturbances  along  shocks  and  Mach  waves.  When  a  generator  Is  new,  the  (rela¬ 
tively  small)  roughness  comes  Inherently  from  the  unevenness  between  the  seg¬ 
mented  electrode  surface  and  the  adjacent  Insulation  (Fig.  1.1);  this  rough¬ 
ness  Is  periodic,  and  extends  over  the  entire  length  of  the  channel.  In  test 
channels  designed  for  the  measurement  of  conductivity  of  a  flowing  plasma.  It 
may  be  Isolated  In  the  middle  of  the  channel  (Fig.  1.2).  After  some  extended 
use,  the  generator  will  show  deeper  and  Irregular  wall  roughness  due  to  erosion, 
corrosion,  burning-out,  and  Inter-electrode  sparks.  As  Indicated  schematically 
In  Figs.  1.1  and  1.2,  shocks  and  Mach  waves  appear  In  a  supersonic  flow,  which 
Interacts  with  the  turbulent  boundary  layer;  and  their  Influence  penetrates  Into 
the  core  flow.  As  a  result,  pressure  rise  In  the  flow  direction  Is  observed  in 
actual  MhJ  generators  with  supersonic  Inlet  conditions  [Ref.  12  and  13].  In  the 
literature,  there  Is  Insufficient  understanding  of  phenomenon  just  described. 

An  Investigation  Is  needed  to  establish,  by  steps  of  Increasing  sophistication, 
an  analytical  model  that  will  provide  additional  understanding,  add  In  the  control 
or  reduction  of  these  effects,  and  contribute  to  the  technology  base  required  for 
development  of  more  efficient  generators. 

As  a  preliminary  endeavor  In  this  direction,  the  present  work  explores  the 
possible  extension  of  Crocco's  theory  of  pseudo-shocks  [Ref.  3]  so  as  (a)  to 
devise  a  simple  analytical  model  for  the  Interaction  near  the  wall  and  Its  Influence 


on  the  core  flow  through  a  supersonic  (high  power-density)  generator,  (b)  to 
explain  quantitatively  (but  on  a  rather  primitive  level)  the  possible  pressure 
ramp  in  the  flow  direction  as  observed  In  actual  generators,  (c)  to  point  out 
the  kind  of  empirical  data  needed  In  completing  or  Improving  the  analysis,  and 
(d)  to  assess  the  role  of  the  model  In  future  studies.  Crocco's  theory  Is  gen¬ 
eralized  to  treat  MHD  pseudo-shocks  as  an  exhibition  of  wave  Interaction;  h1( 
control -volume  analysis  Is  enlarged  to  Incorporate  all  the  MHD  aspects  Including 
load  and  Hall  effects. 

Extensive  numerical  experiments  on  the  solution  of  the  resulting  system  of 
nonlinear  ordinary  differential  equations  demonstrate  actual  qualitative  link 
between  pseudo-shocks  and  pressure  ramps.  Rough  conclusions  are  also  drawn  as 
to  the  various  roles  played  by  the  many  parameters.  The  report  closes  with 
certain  suggestions  for  future  Investigations. 

1.2  CROCCO'S  PSEUDO-SHOCK  THEORY 

In  the  literature,  researchers  In  the  field  over-simplify  the  situation  by 
attempting  to  explain  the  observed  pressure  rise  by  way  of  a  single  normal  shock 
[Ref.  12  and  13]  In  the  core  of  the  generator.  As  a  matter  of  fact,  the  solid 
curves  In  Fig.  1.3  would  be  exactly  the  pressure  variations  In  the  core,  should 
a  single  normal  shock  appear  In  each  case.  The  predicted  variations  are  ob¬ 
viously  too  rough  even  when  viewed  qualitatively.  The  most  outstanding  feature 
of  the  measured  pressure  rise  (Figs.  1.3  and  1.4)  Is  that  the  Increase  toward  a 
maximum  Is  rather  gradual,  reflecting  complicated  Interaction  between  different 
kinds  of  waves  (which  are  oblique  and  hence  comparatively  weak)  before  they 
penetrate  Into  the  core  region.  This  slow  build-up  Is  termed  a  pseudo-shock 
process  In  aon-MHD  gas  dynamics  (Fig.  1.5)  by  Crocco  [Ref.  3].  The  starting 
point  of  Crocco's  work  Is  a  control -volume  analysis  which  recognizes  the  non- 
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uniformity  In  going  from  the  wall  region  to  the  core  (Fig.  1.6).  Some  semi- 
empirical  elements  are  then  elaborated  upon  to  make  the  analysis  self-contained 
[Ref.  8,  14  and  15]. 

For  the  present  Investigation,  the  earlier  and  less  elaborate  (and  hence 
not  self-contained)  phase  of  Crocco's  theory  [Ref.  3]  Is  generalized  to  treat 
MHD  pseudo-shocks.  In  this  approach,  emphasis  is  on  the  qualitative  side. 

Crocco's  control -volume  analysis  is  enlarged  under  the  usual  MHD  approximation 
that  the  magnetic  field  modifies  the  flow,  yet  the  flow  does  not  affect  the 
field.  In  Chapter  2,  mass,  energy  and  momentum- force  balances  are  performed. 

In  Chapter  3,  governing  parameters  are  delineated;  and  in  Chapter  4,  numerical 
examples  are  displayed.  It  Is  to  be  emphasized  that  the  approach  inherently 
Ignores  the  details  of  the  wave  interaction.  From  this  stems  the  simplification 
(and  power)  of  the  analysis.  But,  at  the  same  time  one  learns  to  expect  certain 
lack  of  Information. 

In  the  rest  of  this  section,  we  will  summarize  Crocco's  non-MHD  theory  in 

preparation  for  our  later  generalization. 

* 

Referring  to  Fig.  1.6,  we  have  for  mass  balance 
m  =  m'  +  m"  =  m^ 

and,  for  Invlscld  flow  with  common  pressure  In  the  two  regions  ( boundary- layer- 
and  jet-like), 

^  (p'u'2A'  +  p"u"2A")  =  -(^)A 

and,  finally,  the  energy  balance  (neglecting  wall  heat-transfer  and  viscous 
energy  dissipation) 

2  i2 

^[p'u'A'(h  +  ^  )  +  P"u"A"(h"+^-)]  =  0 

it 

Please  note  that  the  symbols  here,  following  [Ref.  3],  are  not  In  keeping  with 
the  LIST  OF  SYMBOLS. 


The  spatial  coordinate  x  In  the  above  can  actually  be  Ignored  since  the  entire 
system  can  be  rewritten  In  purely  algebraic  form  [Ref.  3]. 

We  may  note  also  that  there  Is  no  entropy  production  In  the  core  region  at 
all;  and  that,  In  the  wall  layer,  the  nonuniformity  across  the  Interface  causes 
entropy  to  be  produced.  The  solution  of  Crocco's  system  Is  a  relatively  simple 
matter;  for  details,  the  reader  Is  referred  to  [Ref.  3],  (Incidentally,  the 
formulation  In  Chapter  2,  of  course,  contains  Crocco's  system  as  a  special  case.) 


CHAPTER  2  FORMULATION  OF  PROBLEM 

2.1  INTRODUCTION 

This  chapter  contains  a  discussion  of  magnetohydrodynamic  channel  flow 
with  a  slowly  varying  cross-section  area.  The  quasi-one-dlmensional  approxi¬ 
mation  will  be  used  to  study  the  phenomenMof  supersonic  flow  under  the 
influence  of  channel  wall  roughness,  In  the  presence  of  Lorentz  force  and 
Joule  heating.  In  order  to  keep  the  discussion  within  reasonable  bounds,  we 
will  consider  only  the  steady  and  Invlscld  flow,  with  no  heat  transfer. 

2.2  GENERAL  RESTRICTIONS 

Before  we  proceed,  let  us  first  limit  the  problem  to  its  acceptable  and 
reasonable  range.  In  order  to  use  the  quasi-one-dlmensional  approach,  the 
cross-section  area  of  the  channel  must  vary  gradually  so  that  the  obliquity  of 
the  wall  with  reference  to  the  axis  Is  always  small.  Therefore,  the  velocity 
component  which  Is  perpendicular  to  the  axis  Is  negligible.  We  also  neglect 
viscous  stresses  and  the  axial  heat  conduction. 

Considering  a  channel  as  shown  In  Fig.  2.1.  The  channel  wall  Is  a  combi¬ 
nation  of  electrical  Insulators  and  conductors.  The  shaded  portions  are  made 
of  electrically  conducting  material,  l.e.,  pairs  of  electrodes  shorted  diagonally. 
The  rest  of  channel  Is  made  of  electrically  insulating  material.  The  mode  of 
operation  Is  the  so-called  "single-load  two-terminal  generation"  (see  Fig.  2.2). 

Further  assumptions  are  Introduced  as  follows: 

(1)  The  temperature  and  pressure  are  approximately  uniform  over  the  channel 
cross-section. 

(2)  There  Is  no  wall  friction. 

(3)  There  Is  no  heat  transfer  at  the  wall. 


2.3  SINGLE-REGION  FLOW 


Let  us  first  consider  the  single-region  MHO  channel  flow  without  bringing 
the  wall  roughness  effect  Into  the  problem. 

Hass  Balance: 

puA  *  *  constant  m  (2.1) 


where  subscript  1  denotes  the  initial  condition,  i.e. ,  the  quantities  measured 
at  the  entrance  of  the  channel. 

Momentum  Balance: 


pudua.d£+  j  B 

dx  dy 

where  j  *  curent  density  in  the  y  direction 

B  *  applied  magnetic  field  strength 

Energy  Balance: 

pjiJH  =  j  e  +  j  f 
Jy  y  Jx  x 

u2 

where  H  *  the  total  enthalpy,  cpT  +  j 


(2.2) 


(2.3) 


j  *  current  density  in  x-dlrectlon,  or  Hall  current 
Ex  *  electrical  field  strength  In  x-dlrectlon,  or  Hall  field 
Ey  ■  electrical  field  strength  in  y-dlrectlon 
Equation  (2.3)  can  also  be  written  in  the  following  form 

*«&T>  &  ♦  &  <T>]  ■  JyEy  ♦  V«  (2-4) 

and  Ex  ■  ctEy  (2.5) 

* 

where  a  »  ratio  of  Hall  field  to  Faraday  field 
and  K  *  ^  (2.6) 

where  K  ■  loading  factor. 

Please  note  that  the  angle  of  Inclination  +  In  Figs.  2.1  and  2.2  is  set  to 
tan’^  a. 


From  the  generalized  Ohm's  law  [Ref.  2],  we  have  two  equations  for  j 


and  J  : 


»  —~7  IX  -  B(E  -  uB)] 
x  1  +  B  x  y 

L  ‘  -S—7  IX  -  uB  +  BE  J 


(2.7) 


(2.8) 


y  1  +  B  y 

where  a  *  electrical  conductivity 
B  *  Hall  parameter 

From  Eqns.  (2.5)  and  (2.6),  Eqns.  (2.7)  and  (2.8)  can  be  rewritten  as 


1  +  B‘ 


Jy  ,  t 
y  1  +  B 


[aK  -  B(K  -  1 )](uB) 


[oBK  +  (K  -  1 )](uB) 


(2.9) 


(2.10) 


Substituting  Eqn.  (2.10)  Into  Eqn.  (2.2),  we  get: 


Hx  *  "  3x  +  +  (K  ~  1 


)](uB‘) 


(2.11) 


Substituting  Eqns.  (2.5),  (2.6),  (2.9)  and  (2.10)  into  Eqn.  (2.4),  we 


PU^Cp  $T+  U  " — g’JCaK  -  B(K  -  l)](uB)'aKuB 


+  ( — - — y)[aBK  +  (K  -  1  )](uB)*KuB 
1  +  B^ 


From  the  Ideal  gas  law. 


and  c 


P  y  - 


(2.12) 


(2.13) 


where  R  -  gas  constant 

y  B  heat  capacity  ratio 

Equations  (2.11)  and  (2.13)  then  reduce  Eqn.  (2.12)  to  the  following: 


*  jf  *  (— *-7>C<.zKz  ♦  (K  -  1)Z3(UBZ)  (2.14) 

’  1  +  6 


2-»/..B24 
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Equations  (2.1),  (2.11)  and  (2.14)  are  the  so-called  mass  balance,  momen¬ 
tum  balance  and  energy  balance  equations  for  single-region  MHD  flows. 

2.4  TWO-REGION  FLOW 

When  a  stream  of  high  temperature  plasma  passes  through  nozzle  and  then 
enters  the  channel  with  high  speed.  It  will  form  a  layer  near  the  wall  because 
of  the  wall  roughness.  In  that  layer  (which  may  be  thick  at  off-design 
conditions),  the  velocity  of  the  plasma  Is  considerably  slower  because  of  the 
wall.  We  call  such  a  layer  the  wall  layer.  We  also  call  the  main  plasma 
flow  In  the  central  portion  of  the  channel  the  core  region.  Thus,  there  are 
two  regions  In  the  MHO  channel.  For  simplicity,  we  add  now  the  following 
assumptions: 

(1)  Wall  regions  near  anode  and  cathode  sides  are  the  same  (or  can  be 
lumped  together). 

(2)  No  pressure  change  occurs  across  a  section  from  region  to  region 
(as  In  unconfined  jets,  or  boundary  layers). 

(3)  There  are  no  friction  and  heat  transfer  between  regions. 

(4)  In  the  regions,  no  property  changes  In  the  z-dlrectlon. 

(5)  In  keeping  with  general  quasl-one-dlmenslonality,  Interfaces  between 

regions  are  supposed  to  be  rather  flat  so  that  Its  normal  everywhere  Is  almost 
* 

In  the  y-dlrectlon. 

Now  we  can  write  down  the  equations  which  govern  the  two- region  flows. 

Mass  Balance: 

PCUCAC  +  ■  P1U1A1  *  constant>  m  (2-15) 

where  subscripts  c  and  ,w  denote  the  core-region  quantities  and  the  wall-layer 
quantities  Individually. 

This  assumption  Is  needed  In  formulating  electrical  conditions  at  the  Inter¬ 
faces. 
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Fig.  2.4  Control -volume  in  the  core. 
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Momentum  Balance  : 

<K<PcucAc  *  PA>  -  -  £  ■ *0  -  $  \  +  V  *  W  (2j6) 
where  A  *  cross-section  area  of  core-region 
■  cross-section  area  of  wall-layer 

Let  A  be  the  total  cross-section  area  of  the  channel  flow.  Then, 

A  *  Ac  +  Aw.  Together  with  Eq.  (2.10),  Eq.  (2.16)  becomes: 


Before  we  can  go  further  to  study  this  two-region  problem,  we  have  to  find 
the  key  point  which  connects  the  core  region  and  the  wall  layer.  To  achieve  this 
goal,  let  us  first  examine  the  core  region  energy  equation. 


where  i(p  u  A  )  5  0  for  the  core  region  Is  to  be  In  force  throughout  the  present 


Investigation. 

•far  later  use,  we  will  also  note  here  that  the  mass-rate  out  of  the  core  across  Its 
lateral  surface  between  x  and  x+dx  Is  pcucAc  -  CpcucAc  +  ^PcucAc)‘dx^  s 

-  ^PcucAc)*dx  •  that  1*»  -  ^pcucAc)  per  unit  axial  length  (see  Fig.  2.4). 


It  Is  Important  to  realize  the  full  Implication  of  the  above  restriction, 

-r-(p  u  A  )  *  0.  Under  this  restriction,  the  fluid  particles  in  the  core  region 
ax  c  c  c 

will  keep  on  carrying  energy  out  of  It,  and  entering  the  wall  layer.  If  ^ecucAc) 
were  positive,  energy  would  be  carried  into  the  core  region  from  the  wall  layer. 

In  that  case,  the  contribution  term  (the  second  |erm  on  the  lefthand  sid^)  of 

Eq.  (2.19)  would  have  to  be  replaced  by  -{h^+  T^pcucV*  or  (hw  +  ’F^pwuwV' 

This  can  be  expressed  by  the  following  statement: - 

(Net  energy  carried  out  across  sections  at  x  and  x+dx)  +  (Energy  carried  out  of  the 
core  across  the  lateral  surface)  ■  (Electromagnetic  work  done  on  the  plasma  In  the 
core  between  x  and  x+dx.) 

Expanding  Eq.  (2.19),  together  with  Eqs.  (2.9)  and  (2.10),  yields 
2  2  2 
<hc  +  *  VAsA  *  T>  •  (hc  +  t’&VcV 

-  j(^f)[.eKc- sc(Kc-  1)](UCB)-(.CKC«CB)» (^r)C«e»ei‘c+ <v 
c  c 

■  (it7)Uckc  va  Kc  ■ 

c 

l.e..  ♦(7^7)[«cltc+C'c-,>2^c!|2>  '2-20> 


where  Eqs.  (2.11)  and  (2.13)  have  been  employed. 

Similarly,  we  can  derive  the  momentum  equation  of  core  region  as  follows: 

lr(pcucAc)'uc^pcucAc)  *  "  ^Ac+(7^7)[acBcKc+(Kc"1)](ucB)2,Ac 


Expanding  the  above  equation,  we  get 


“c^^cAc)+  (pcucAc )_ar  •  ‘JcJrK“cAc>  =  ■  H&  + 

‘77T?>t“cBcKc  -  ‘KC  -  i)3(UcB2)Ac 


or,  using  Eq.  (2.13), 


“  -  ^+(7^7)[acecKc+(Kc-1):,(ucB2)  (2‘21) 


Comparing  Eqs.  (2.21)  and  (2.20)  to  Eqs.  (2.11)  and  (2.14),  it  is 

obvious  that  the  core-region  formulas  are  similar  to  the  single-region  ones.  The 

wall-layer  equations  are  not  similar  since  the  quantities  carried  into  the  wall- 
2  2 

ur  IT! 

layer  are  hc  +  -Tp-  and  uc,  not  hw  +  and  i^. 

For  a  two-region  flow,  the  core  region  cannot  be  studied  in  a  decoupled  manner 
(Independently  of  the  wall  layer)  since  o  and  K  are  related  to  wall-layer  quanti- 

C  C 

ties  through  electrical  conditions  at  the  Interface.  We  will  make  further  dis¬ 
cussion  on  this  point  later  on. 

Now  that  we  have  the  equations  governing  the  two  regions  comblnedly  and  the 
core-region  equations,  we  can  then  easily  obtain  the  wall-layer  equations  by  sub¬ 
tracting.  For  example,  subtracting  Eq.  (2.20)  from  Eq.  (2.18)  will  yield  the 
energy  equation  of  the  wall  layer: 

2  2 

&pwuwAw(hw  +  T)]  +  (hc  +  T}  ^(pcucAc)  =  (jxwExw  + jywEyw)Aw  (2*22) 

To  continue  our  analysis,  we  need  the  assumed  division  of  masses  in  the  two 
regions  In  the  spirit  of  Crocco  [3],  Let  us  introduce  the  mass  ratio  y  such  that 


■  uGim 

■>cucAc  ■  0  - 


(2.23) 


where  m  »  total  mass  flow  rate 


, .  _  mass  flow  rate  In  wall  layer 
W0Q  total  mass  flow  rate 

d  dgM 

Following  our  previous  restriction,  ^fcucAc)  -  0»  we  see  that  -3"-  t  0. 

An  alternative  approach  would  be  to  start  with  a  division  of  flow  cross-sectional 
areas  In  the  two  regions: 

*w  ■  v&va 

AC  ■  0  -  '■ft)3'A 

"here 

It  may  appear  superficially  that  \>fy  should  be  more  directly  measurable. 

But,  to  guarantee  ^  >  0,  It  Is  safer  to  start  with  vfy  given.  And  can  then 
be  calculated  as  a  part  of  the  result. 

In  eq.  (2.16),  we  have  already  employed  the  relation  aQQ  *  Aw  +  Ac •  Together 
with  Eq.  (2.23),  this  relation  now  yields 

“to  1 1  *  “06 .  a6) 

“A  “c“c  ■ 

Introducing  further  the  Ideal  gas  law,  we  have 


(2.24) 


From  Eq.  (2.22),  we  can  write 

dh. 


.  dh  du  11  a 

hw  ^r(pwuwAw)  +  <pw\Aw>  ~3x"+  <pwuwAw)uw  dx  +T"  dx^pwuwAw^ 


+  ^hc  +  7^  33T(m-pwuwAw)  *  ^xwExw+ ^ywEyw^ 
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2 

dh  (u  -  yJ  ,  a  .  .  .  .  dp  . 

^pwmwAW^"3x*+  ^hw"  hc - 2  ^3x  <pwuwV  *  “w**  dx 


<7^>^*V1>XB‘Aw 


2-,  2B2, 


1  +  e 


that 


is. 


Finally,  replacing  p^u^  by  w(x)-m  yields 


Equation  (2.26)  can  also  be  rewritten  In  terms  of  v(x) : 


(2.27) 


Equations  (2.20),  (2.21),  (2.24),  (2.27)  and  (2.28)  then  govern  the  MHO  two 
region  flow.  Before  we  proceed,  we  must  note  carefully  the  dependent  variables 
In  the  above  five  governing  equations.  Let  us  list  all  the  dependent  variables 

here: 

p  »  p(x) ,  uc  •  uc(x),  •  uw(x),  Tc  ■  Tc(x),  Tw  -  Tw(x). 

The  auxiliary  quantities, 
oc  ■  fct.  (Tc,p) 
ow  ■  fct.  (T^,p) 


3.4 


$c  *  fct.  (T.,p,B) 

0*  -  fct.  (Tw,p,B) 

are  seen  to  be  related  to  the  other  quantities  via  known  functional  forms  (to 
be  quoted  later). 

The  given  quantities  of  the  problem  are  A(x),  B(x),  aw,  1^,  y(x),  and  the 

Initial  values  of  p,  u  ,  u  ,  T  and  T  .  (a  and  K  will  be  related  to  a  and 

w  c  w  c  c  c  w 

In  Section  2.5.)  In  principle,  aw  =  tan  i  (see  Fig.  2.1)  can  be  a  general 
function  of  x;  although  a  constant  will  be  used  throughout  the  present  study. 

It  must  be  emphasized  that  y  as  a  given  function  of  x  reflects  very  clearly 
the  qualitative  nature  of  the  present  theory.  The  formulation  yields  a  solution 
for  a  prescribed  but  In  a  truly  self-contained  framework,  y(x)  ought  to 

come  out  as  a  part  of  the  solution.  As  It  stands,  one  will  have  to  speculate 
about  the  functional  form  of  y(x)  and  calculate  examples  that  are  hopefully  repre¬ 
sentative  of  actual  cases  in  a  qualitative  manner.  It  Is  also  not  truly  satisfying 
to  view  the  theory  as  being  semi-empirical  (In  the  sense  that,  for  every  actually 
measured  y(x),  one  can  predict  a  complete  flow  situation),  simply  because  (1)  It 
Is  Impossible  to  measure  locally  the  relative  mass  flow  In  the  two  regions,  and 

(2)  the  quasi -one-dimensional  approximation  itself  obliterates  the  true  meaning 
* 

of  an  Interface  .  Of  course,  such  difficulties  are  already  an  Integral  part  of 
Crocco's  classical  pseudo-shock  theory  [ref.  3];  we  have  only  Inherited  Crocco's 
qualitative  legacy  In  our  attempted  generalization. 

2.5  ENTROPY  PRODUCTION 

Applying  the  perfect  gas  law  and  the  general  relation  between  entropy  .8, 
Internal  energy  U  and  volume  V,  one  has 


This  Is  quite  similar  to  the  fact  that  in  the  boundary  layer  approximation  the 
"edge"  of  the  boundary  layer  Is  a  hazy  concept. 


J 


jo  _  dU  .  l>dV 

dS  -  T+  “’T 


where  cp  *  heat  capacity  at  constant  pressure  (*  ^^-) 


Therefore, 


or, 


(2.29) 


1  -  (  1  \  1  dT  1  dt*  ,,  ,nA 

*35"  (7^T}  T3x  "  F  dx  (2*30) 

By  comparing  Eq.  (2.30)  with  Eqs.  (2.20)  and  (2.27),  we  may  have  the  equations 

of  entropy  change  In  the  core-region  and  the  wall  layer  Individually: 

1  <»c 


1-3T  “  ^7~?^acKc+ (Kc‘ 1)2](ucb2) 

1  ec 

ii?  •  ^(777)Cowkw+  (kw*1)2](uwb2) 

Bw 

(u  37)(-RT")C(A)(Tw'Tc)  ’  ^uc'uw)2] 


(2.31) 


(2.32) 


where: subscript  c  denotes  the  quantities  In  the  core-region, 
subscript  w  denotes  the  quantities  In  the  wall-layer. 

It  Is  thus  seen  that  the  entropy  production  In  the  core  Is  due  solely  to  the  Joule 
heating  (which  Is  a  form  of  energy  dissipation).  In  the  wall-layer,  there  are  two 
sources  of  entropy  production:  Joule  heating  and  non-uniform  distributions  of  u 
and  T  In  the  two  regions. 

For  the  non-magnetlc  case  [Ref.  3],  the  core  flow  Is  Isentroplc  and  the 


entropy  production  In  the  wall-layer  Is  due  only  to  the  non-uniform  distributions 
In  the  two  regions. 


2.6  ELECTRICAL  CONDITIONS  AT  THE  INTERFACE 


In  the  previous  discussion*  we  encountered  several  electrical  quantities: 

J  ,  j  ,  E  and  E  .  To  find  the  relationships  which  can  link  these  electrical 

Jr  X  jr  X 

quantities  all  together  Is  the  key  to  the  solution  of  our  problem.  According  to 
Sec.  1.6  of  [Ref.  4],  we  have  (Fig.  2.5): 

<n  x  ■  0 

<n  •  *  0 


where  <  >  stands  for  jumps  across  the  Interface: 


E  »  E 
xc  xw 

(2.33) 

j  ■  j 

Jyc  Jyw 

(2.34) 

(Note  that  f  end  J„c  !•]„.) 

By  using  the  Interrelationships  expressed  In  Eqs. 

(2.5),  (2.6),  (2.34) 

and  (2.35),  one  has 

“w  _  Exw/Evw  _  Evc 
*c  £ 

(2.35) 

"w  _  Aw'c, 

r  tVfuff)  'a  ,(u  ' 
c  yc  'c'  w  “w 

(2.36) 

or,  E  u  u 

\  ■  w  ■ 

(2.37) 

where  K*  ^  Is  a  given  quantity.* 
c 

Combining  Eqs.  (2.10),  (2.36)  and  (2.34)  will  yield 


(^)[VA+(Vl)](UcB)  *  ("\)[VwV(Kk'1)](uwB) 
1  +  6-  1  +  6*. 


In  usual  designs,  )(■  0.5  which  Is  an  optimizing  value  In  some  sense. 


which  can  also  be  written  as  follows: 


( K  u  s  (_!w^)(K  1)u  (_!c^)(K  1)u 
i  +  e:  i  +  C  i  +  e:  w  w  i  +  e:  c  c 


(2.38) 


We  also  mention  here  as  an  addendum  that,  from  the  literature 


B  «  B/T/P 


(2.39) 


and  o  «  Tm/Pn  (2.40) 

where  m  and  n  are  universal  constants  (depending  on  the  temperature  range) 
for  all  plasmas. 


CHAPTER  3  DIMENSIONLESS  PARAMETERS 


3.1  NONDIMENSIONALIZATION 

Applying  the  process  of  dimensional  analysis  to  the  five  governing  equations, 
Eqs.  (2.20),  (2.21),  (2.24),  (2.27)  and  (2.28),  we  Introduce  the  following  repre¬ 
sentative  quantities: 

Density  —  f>c1 
Velocity  —  ucj 
Area  —  A^ 

Generator  Channel  Width  —  b^,  see  Fig.  2.1 
Magnetic  Field  Strength  — 

Pressure  —  ^ 

Temperature  —  T^ 

Conductivity  —  o  ^ 

where  subscript  1  denotes  the  Initial  conditions,  l.e.  the  quantities  at  the 
entrance  of  MHD  channel. 

Then,  define  the  following  dimensionless  variables: 


where  prime  denotes  dimensionless  quantities. 

Substituting  the  above  dimensionless  quantities  into  the  governing  equa¬ 
tions  derived  In  Chapter  2,  we  then  have  the  new  governing  equations  with  dimension 
less  forms  as  follows: 


(1)  From  Eq.  (2.20), 


t>>  v  d(T.T')  d(t>>)  ariat'  9  9  5  19 

<t^)(7^t)  d(b;x'j  “  Ttsprr f  (777)C“cKc* (Kc‘ 11  ](uciuJ8i 


2B'2) 


After  simplifying,  we  have 


<r><7?T>  ST ■  3F +  (->rcA»7^)[^Kc* <V (3.1) 


(2)  From  Eq.  (2.21), 

f*>u  .u*  d(u  .u' )  d(l>>)  o.o l 

(  5T  T1 — )  ~a1Z"Z'\ - m-?'V  +  (■ - B^K  + 


n?T^Tp  _3^T  =  ‘  dTFjx’T  *  (7^!)CVckc+  (KC‘  1)](uclucBiB'2) 

uL  du'  Hb,  b4°r<UMB4  ar  9 

-a*  -  -  +  (-^UU)(C  }[  K  (K  1)](u,b,2) 


KV  TT  3x 


P1  )7f  CCC  v,Nc 


Observing  the  above  equation,  let  us  multiply  a  quantity  y  (the  ratio  of  heat 
capac11&)on  both  denominator  and  numerator  In  lefthand  side  and  notice  that 


is  the  Mach  Number  at  the  Inlet  of  the  channel  In  the  core  region. 


Then,  we  have 


a#'  -  ar ♦  !  <3-2> 

(3)  From  Eq.  (2.27), 

f,V'  w  Y  \  d(TClTW}  /  V'  vrfyR  xT  (T,  T,x  UCl/u,  ,  x2,  dy 

^  d7^j7T+  (RT^n;)C(:rT)  Tc1(Tw‘Tc)‘“nuc  V  %pr*) 


)  Cf  iO  0  0  0  0  9 

■  “HH^T +  V  d2](Uc1U;b2b,: 2' 


0 

>'(Tj,(7^T)  3^’  4j)[(7^T,(Ti‘  Ti> '  (-ri)(“;-u;)2l|^ 

1  1  Bw 

(4)  From  Eq.  (2.28), 

|>  t>*  d(u  u* }  {>  )>' 

p(irr^rj)(Uctui)  d(Bc;xrV+  (^T;)(uci^)(ui-  uP*nyn- 

- “3is{xT+  “<777>tv,v <s.- ’)](“ci“iBfB' ) 

«  fr'u'  du‘  «  f>V  . 

u{ylii){lf)  d3T+  UP  ^  3 

-■'a?-*  u(  -  cp  cf  f')(777)CVwKw+  (V1)](uiB‘2) 

*  pw 

(5)  From  Eq.  (2.24), 

.  u  ,u'  u  .u'  u  . u*  u  .  u'  u' 


(3.3) 


(3.4) 


*  -^+^‘J)]  (3.5) 

C  w  C  « 

^ElUc1®1 

Here  let  us  Introduce  a  new  parameter,  S,  ■  —  h--  ,  called  the  Interaction 

l  r1  - 

Parameter.  To  summarize,  now  we  have 

(Tj»7^r>  a^‘  hf+  si(TrV,':'cl'ct(Kc-1)Jl("i8'2)  (3-6) 

0 

9  Ni  du’  AI  o’  « 

{T'ii){’TT)W  "■axr+S1(777)CacecKc+  (Kc"  1)3(uiB'  ) 


(3.7) 


3* 


“3T  *  “si<^~?)Cow,vt  0^-  >)2](u;b,2) 

-  K'  du'  «  ^'u'  , 

M(yMct)(-1y-)aF-+  (Y^,)(-yL)(,t-  uc^T  = 

w  w 

•  0* 

-  -ar ♦  “Vt^AVV d](u;b'2) 
>'*  'rwilW-iW-Wi 


(3.8) 


(3.9) 


(3.10) 


If  we  take  the  natural  logarithm  on  both  sides  of  Eq.  (3.10),  It  will  yield 
*"fr'  «  tn[T'u;  +  w(ru' -  TV)]  -  in(A'u^) 

Then,  differentiating  the  above  equation  with  respect  to  x',  we  have 

l  ,4>'  f  du'  dT*  du'  dT'  du'  dT'  / 

r  ar  a{C(Tc  d^ +  u;  ^ +  ^  "c  w-  t;  u; 

A“i +  “(TW-wf  (3.11) 

Eqs.  (3.6)  through  (3.9)  together  with  Eq.  (3.11)  form  a  system  of  non-linear 
differential  equations,  which  can  be  solved  numerically. 

3.2  DIMENSIONLESS  FORMS  OF  AUXILIARY  EQUATIONS  AND  CONDITIONS 

Referring  to  Sec.  2.6  and  Introducing  the  dimensionless  quantities  defined 
In  previous  section,  we  can  get  the  following  relationships  which  connect  the 
loading  factor  K,  the  velocity  of  plasma  u  and  a  (the  ratio  of  Hall  to  Faraday 
^leld): 


"  Klif> 


(3.12) 


'  Vtyty  ■  «£> 


I 
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or-  »e  ■  %<£> 

Substituting  Eq.  (3.13)  Into  Eq.  (2.39)  and  nondlmenslonallzlng  will  yield 


(3.13) 


O  (j*  ft 

/  Cl°CPC  Cl  W  Ww  *  \/  w  ..  -i\ 

( - 5 - t-k?“)(<*wK-u  4u') 

1+4  1 + s,  c  1 


ac1ffw6w. 


In  Eq.  (3.14),  there  are  two  quantities  which  vary  with  respect  to  x',  l.e.,  B  and 
o'.  As  mentioned  In  Sec.  2.6,  Is  related  to  the  mth  power  of  temperature  and 
the  nth  power  of  pressure.  And  from  the  literature,  the  values  of  m  and  n  are 
respectively  10  and  0.5.  Introducing  these  values  to  Eq.  (2.41)  and  rewriting 
Eqs.  (2.40)  and  (2.41)  In  dimensionless  forms,  one  has 


°c1  'cl  /  P1 


In  later  calculations,  we  will  let  |!J^  equal  to  1 
become  the  following  forms: 


(3.15) 

(3.16) 

Then,  the  above  equations 


B'VT 


(3.17) 


(Tl)10 

Jr 


(3.18) 


In  existing  diagonal-wall  generators,  Bci  Is  always  close  to  1;  there  is  also 
rational  ground  for  scj  -  1  In  good  designs,  see  [Ref.  7j. 


We  also  change  the  electrical  quantities  to  dimensionless  forms  by  the  following 


steps: 


Jyw '  777  ■VA* ( 


Similarly, 


Jxc  ■  77?  CacV  6c(l!c',)](uclW) 

'  ^°ciuc1B1^7~”!?^0cKc~  ffc^Ec"  E^ucB*  ^ 
1  Bc 

^xw- 


(3.19) 


(3.20) 


(3.21) 


From  Eqs.  (2.6)  and  (2.5),  let  us  rewrite  the  electrical  field  strength  as  follows: 
In  the  y  direction,  we  have 


Eyc  *  KcuclW  ■  <“c1B1><W> 

Eyw-<“c,V<W> 


In  the  x  direction,  E  «  E 

w  AW 


Exc  *  °cEyc  ■  <uC1B1>**cW> 

Now  let  us  define  the  following  dimensionless  electrical  quantities: 


(3.22) 

(3.23) 


(3.24) 


41  «  41  «  ^XC.  {I  .  ^XW 

°cl“clV  xc  »e1uciV  *  °c1uc1B1* 


E^-^  Eic-^ 


And  from  Eqs.  (2.34)  and  (2.35),  ■  J^c,  EJC  -  E^.  From  above  definition*, 

Eqs.  ^ 3. 1 9 )  through  (3.24)  become 


v 


3^ 


(3.25) 

(3.26) 

(rrSH V»w<Kw-,>«‘CB'> 

1 

(3.27) 

e;c  ■  W 

(3.28) 

E*  -  w  - 

(3.29) 

E*<*  E»>  ‘“Ac 

(3.30) 

With  these  dimensionless  electrical  quantities  on  hand  ,  we  can  find  the  trend  of 
these  variables  In  a  two-region  MHD  channel  flow  problem. 

Another  object  which  we  are  interested  In  is  the  area  variation.  In  order 
to  develop  the  area  ratio,  we  combine  Eq.  (2.15)  and  Eq.  (2.23)  and  get 

-.“A '  “VciN 


Simply  by  mathematical  manipulation,  we  have 


*; 


(3.31) 


where  ■  the  ratio  of  wall-layer  cross-section  area  to  the  cross-section  area 
at  entrance,  l.e.,  A^/A^ 


Similarly,  T' 

A£  »  (1  -  M)(^jrpr-) 


(3.32) 

or,  AJ  -  A*  -  A;  (3.33) 

From  Eqs.  (3.31)  and  (3.32),  we  see  that  the  cross-section  area  of  wall -layer  Is 
proportional  to  temperature  and  mass  ratio,  and  inverse  to  velocity  and  pressure. 


3.3  PHYSICAL  SIGNIFICANCE  OF  GOVERNING  PARAMETERS  Mci  AND  S< 

In  Sec.  3.1,  there  are  two  parameters  which  appeare  In  governing  equations, 
Eqs.  (3.6)  through  (3.9),  that  play  very  Important  roles  In  the  two-region  WD 
channel  flow.  These  two  parameters  are  the  Mach  Number  at  the  Inlet  (Mc1)  and 
the  Interaction  Parameter  (S^).  Since  Mc^  and  S^  almost  appeare  In  each  equa¬ 
tion,  we  can  predict  that  these  two  governing  parameters  w<11  Influence  strongly 
the  behavior  of  plasma  In  the  generator.  For  example.  In  Lq.  (3.7),  we  have  In 
the  x  direction  and  on  an  elementary  plasma  cell: 

(Inertial  Force)  ■  (Pressure  Force)  +  (Lorentz  Force) 

Thus, 

c 

Lorentz  Force  a  1 
Inertial  Force  12“ 

"cl 

Lorentz  Force  a  s 
Pressure  Force  1 

which  means  that  the  larger  the  Interaction  parameter  S^,  the  larger  the  Lorentz 
Force  while  the  smaller  the  Inertial  Force  and  the  Pressure  Force.  On  the  other 
hand,  the  larger  the  Mach  number,  the  smaller  the  Lorentz  Force  while  the  larger 
the  Inertial  Force. 

Also,  In  Eq.  (3.6), 

(Change  In  Internal  Energy)  *  (Work  Done)  +  (Joule  Heating) 

l.e. , 

Joule  Heating  B  s 
Energy  Change  “  1 


Joule  Heating  B  s 
Work  3l 

which  means  that  large  value  of  Sj  will  produce  strong  effect  of  Joule  Heating 
and  vice  versa. 
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For  this  problem,  since  plasma  which  comes  out  from  nozzle  enters  the 
channel  with  supersonic  speed,  the  Mach  Number  M^  Is  always  greater  than  1. 

To  estimate  representative  value  of  Sj,  let  us  use  some  ranges  of  data 
gathered  from  existing  literature: 

—  25  to  60  mho/m 

—  1  to  3  atm 

—  around  0.15  m 
B.j  —  2  to  6  Tesla 

u  i  —  around  750  m/s ec 
T.  —  3000  to  4000  °K 

Therefore,  Is  approximately  In  the  range  from  0.3  to  2.0.  As  a  result  we 
will  use  Sj  *  0.3  ^  0.5  In  our  calculations,  to  keep  to  the  lower  (and  more 
realistic)  end  of  the  range. 
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CHAPTER  4  RESULTS  AND  DISCUSSION 
4.1  NUMERICAL  EXAMPLES 

Several  examples*  are  worked  for  a  plasma  w1th?T>1.1,  to  display  certain 
trends.  Figure  4.2  to  4.14  show  plots  of  various  quantities  for  the  case  where 


A'  • 

1  +  0.2x' 

v  ■ 

0.0003  +  0. 1 2x 1 ,  v  <  0.3 

a 

0.3,  v  *  0.3 

"if- 

1,  u^  -  0.008 

Tc(  ■ 

Til  "  1 

Pi  ■ 

1 

V 

0.4 

K  ■ 

0.5 

a  1 

w 

-1 

The  curves  marked  with  crosses  are  for  Mc^  *  1.4,  and  those  with  small  squares, 

"ci  -  1.5. 

For  comparison.  Fig.  4.1  shows  the  pressure  variation  for  the  single-region 
model  where  A',  Sj,  K,  pj  are  kept  the  same;  o  *  -1,  uj  *  1,  TJ  *  1;  and  *  1.4 
(1.5)  for  the  curve  marked  with  crosses  (squares).  It  Is  seen  that  the  plasma 
pressure  decreases  monotonously  In  the  flow  direction.  This,  no  doubt,  would 
be  Indicative  of  the  case  where  there  Is  no  appreciable  wall  roughness.  In  con¬ 
trast,  Fig.  4.2  shows  that  the  plasma  rises  In  Its  pressure  rather  gradually  In 
the  flow  direction,  exhibiting  a  (WD)  pseudo-shock  when  a  two-region  model  Is 
adopted.  This  lends  credence  to  the  qualitative  rationale  behind  an  explanation 
of  the  observed  pressure  ramp  by  way  of  a  two-region  pseudo-shock.  This  highlight 


All  with  B'(x)  »  1  ♦  0.5x' ,  for  x’  i  1;  and  1.5  (1  -  0.1x‘),  for  x'  >  1. 


of  the  present  Investigation*  however*  does  not  shine  too  much  quantitatively* 
since  the  ratio  of  the  mass  flow  In  the  wall  layer,  u(x'),  has  to  be  conjectured 
for  the  calculation.  Furthermore*  there  does  not  seem  to  be  any  practical  means 
that  could  yield  actual  data  leading  to  n(x')  for  a  given  run  of  a  real  generator. 
Thus,  It  would  be  futile  to  hope  for  a  close  comparison  between  a  calculation 
and  data  recorded  for  an  actual  run;  yet.  It  Is  perhaps  possible  to  gain  some 
understanding  of  the  underlying  mechanism,  and  the  general  trend  by  examining 
the  numerical  plots. 

Figure  4.3  shows  the  local  Mach  number  In  the  core;  and  Fig.  4.4*  that  In 
the  wall  layer.  Mc  Is  seen  to  decrease  monotonously;  first  Increases  and 
then  decreases.  Sharing  the  same  trend  Is  the  velocity  In  the  core  (Fig.  4.5); 
and  that  In  the  wall  layer  (Fig.  4.6).  Thus,  as  the  wall  layer  grows  Its  velocity 
Is  first  promoted  by  the  core  flow  (at  the  expense  of  the  core  velocity);  but 
later  decelerates,  probably  due  to  the  Lorentz  force  accompanying  the  generation 
of  electricity.  The  core  temperature,  as  shown  In  Fig.  4.7,  Increases  monoton¬ 
ously  (but  only  by  about  10*).  The  wall-layer  temperature  (Fig.  4.8)  takes  a 
rather  sudden  jump  near  the  entrance,  and  then  Increases  more  gradually  until 
15  or  30*  higher  than  the  Initial  value.  The  reason  for  the  jump  In  wall-layer 
temperature  at  the  entrance  Is  not  clear;  but  an  educated  guess  Is  that  It  Is 
brought  on  by  the  specific  y(x')  we  used  (either  Its  Initial  value  or  Initial 
slope).  It  probably  has  no  realistic  counterpart  In  actual  cases. 

Figures  4.9  to  4.14  show  various  electrical  quantities.  Since  the  situation 
we  are  Investigating  Is  patently  off  the  design  condition,  we  will  refrain  from 
commenting  on  these  plots  as  they  serve  no  practical  purpose,  except  to  be  complete 
and  to  display  rather  wild  off-design  variations. 
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Similarly,  Figs.  4.15  to  4.23  show  plots  for  the  case  with 
A*  -  1  ♦  0.061x' 

Si  -  0.3 

and  Mc1  *  1.5  (curves  marked  with  triangles),  1.6  (circles),  1.7  (crosses);  all 
the  other  parameters  are  the  same  as  In  the  previous  case.  These  figures  show 
the  same  trend  as  before. 

Another  set  of  results  Is  given  In  Figs.  4.24-4.35.  Here,  =  0.4,  *  1.7; 

everything  else  Is  the  same  as  for  the  previous  case.  In  addition.  Fig.  4.36  also 
shows  the  relative  cross-sectional  area  of  the  wall  layer  (as  fraction  of  the 
total  area)  for  this  case.  It  Is  seen  that  the  wall  layer  thickens  quickly  near 
the  entrance;  and  eventually  fills  more  than  80%  of  the  duct  for  this  particular 
case. 

Finally,  In  Figs.  4.36-4.38,  Influences  on  the  pressure  variation, 
due  to  changes  In  Inlet  Mach  number,  maximum  u-value,  and  degree  of  duct  diver¬ 
gence  are  (respectively)  plotted.  In  all  these,  S1  *  0.4;  and  u  j,  uw1,  etc., 
are  all  the  same  as  before. 

In  Fig.  4.36, 

A'  -  1  +  0. 061  x' 
w  *  0.0003  +  0.1 2x' 

*  0.3  when  the  calculated  u  -  0.3 

The  curve  marked  with  triangles  Is  the  case  with  Mci  *  1.5;  that  with  circles, 
1.6;  and  that  with  crosses,  1.7.  Thus,  the  "strength"  of  the  pseudo-shock  In¬ 
creases  with  Mc1. 

In  Fig.  4.37,  M  ^  a  1.5,  A'(x')  Is  the  same  as  for  Fig.  4.36. 
u  *  0.0003  +  O.lx' ,  x'  <  1 


0.1003  +  0.12(x'  -  1),  x'  >  1 
Umax  when  the  calci|,,ated  v  -  w^x 
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The  curve  marked  with  triangles  represents  the  case  with  u|WX  a  0.2;  that  with 

crossed  circles,  0.3;  that  with  squares,  0.4;  and  that  with  crosses,  0.5.  The 

trend  of  larger  u  Is  seen  to  make  the  pressure  ramp  flatter. 

Figure  4.38  shares  with  Fig.  4.36  the  same  v(x'h  and  with  Fig.  4.37,  the 

same  Mcj.  For  the  curve  marked  with  triangles,  we  have: 

A'  *  1  +  0. 061 x ' 
that  with  crossed  circles: 

A'  -  1  +  O.lx' 
that  with  squares: 

A'  *  1  +  0.2x‘ 
and  that  with  crosses, 

A'  =  1  +  0.3x' 

There  Is  no  clear  trend  discernible.  However,  it  Is  noted  that  the  pressure  vari¬ 
ation  Is  very  sensitive  to  the  duct  divergence. 

4.2  CONCLUSIONS 

Under  the  present  Grant,  we  have  successfully  demonstrated  the  feasibility 
of  extending  Crocco's  pseudo-shock  model  to  explain  qualitatively  the  observed 
pressure  rise  In  a  rather  extended  portion  of  supersonic  MHO  generators,  with 
rough  walls  and  at  off-design  conditions. 

For  future  refinement,  we  may  list  a  number  of  ways  by  which  (singly  or 
comblnedly)  the  formulation  may  be  Improved: 

(1)  Heat  transfer  and  friction  at  the  wall  can  be  introduced  through  addi¬ 
tional  coefficients. 

(2)  Heat  transfer  and  friction  at  the  Interface  can  be  accounted  for. 

(3)  The  formulation  can  be  coupled  with  a  boundary  layer  analysis  (using 
Karman-Pohl hausen  technique)  at  the  wall. 


(4)  The  wall  layers  at  the  cathode  and  anode  can  be  distinguished,  and  the 
formulation  extended  to  a  three-region  model. 

(5)  Cases  with  du/dx  <  0  can  be  tested  (with  due  modification  of  the 
governing  equations). 

However,  the  major  drawback  In  our  analysis  (also  in  Crocco's),  l.e.,  the 
empirical  assignment  of  a  y(x)  to  start  the  calculation  will  always  be  with  us. 

It  may  also  be  Interesting  to  apply  the  two-region  model  to  flow  apparata 
other  than  power  generators. 
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APPENDIX  —  COMPUTER  PROGRAM 


1 


IV  G  LEVEL  21  MAIN  DATE  =  80189  14/30/26 

C  THIS  PROGROM  IS  TO  CALCULATE  THE  DEPENDENT  VARIABLES  OF  TWO-REGION  BHD 
C  CHANNEL  FLOW  .  TWO  SUBROUTINES  ARE  USED  ,  NAMELY  FCT  AND  CAC  INDIVIDU- 
C  ALLY  .  THE  FORMER  USES  THE  MODIFIED  GAUSS-JORDGH  METHOD  TO  SOLVE  FIVE 
C  DIFFERENTIAL  EQUATIONS  TO  GET  THE  GRADIENT  OF  DEPENDENT  VARIABLES  . 

C  THE  LATTER  DEFINES  ALL  AUXILLIARY  QUANTITIES  AND  ELEMENTS  OF  THE  MATR- 
C  IX  WHICH  IS  GOING  TO  BE  USED  IN  SUBROUTINE  FCT  . 

c  the  Parameters  being  used  in  this  program  are  defined  as  follows  : 


C  nil - VELOCITY  OF  CORE-REGION 

C  Y  (2 ) - VELOCITY  OF  WALL  LAYER  REGION 

C  Y  (3 ) - TEMPERATURE  OF  CORE-REGION 

C  Y  (4) -  —  -TEMPERATURE  RATIO  OF  WALL  LAYER  REGION 
C  Y  (5  ) - PRESSURE  RATIO 

C  X - DISTANCE  RATIO  ALONG  THE  X-DIFECTICN 

C  AM - MACH  NUMBER  AT  ENTRANCE 

C  CK - LOADING  FACTOR 

C  s - INTERACTION  PARAMETER 

C  AR - RATIO  OF  HALL  FIELD  TO  FARADAY  FIELD 

C  XP - INTERVAL  SIZE  FOR  PRINTING  OUT 

C  OER Y  (I) — — DERI VATIVESOF  DEPENDENT  VARIABLES 
C  GC - RATIO  OF  SPECIFIC  HEAT  CAPACITIES 


IMPLICIT  REAL*8  (A-H,P-Z) 

DIMENSION  Y  (5» , DERI  (5) 

COMMON  GC,AR,CK, AM,S 
C  DEF  INS  CONSTANTS  : 

GC= 1. 10DU0 
AR  =  -  1. 00D00 
CK=0.  50D00 
AM= 1 . 5UD0  0 

C  DEFINE  INITIAL  VALUES  OF  VARIABLES  : 

Y  (1)  =1.  00D00 

Y  (2) =0. 80D-2 

Y  (3)  =  1 . 00D00 

Y  (4)  =1.  OODOti 

Y  (5)  =1. 00E00 
S=J.40D00/1 .5GC00 
X=0. 00D00 

WRITE (6,300) 

300  FORMAT  ('  1 '  ,  12X,  '  X»  ,  19 X,  '  Y  { 1 )  1  , 16X  ,  •  Y  (2)  •  ,  1 0X  ,  '  Y  (3)  '  r  IfcX,  '  Y  (4)  •  ,16X 
1  ,'Y(5)') 

XP=0 .  OODOO 
H=0. 50D-2 

10  CALL  FCT  ( X , Y , DEL Y) 

IF  (DABS  (X-XP)  -0.  ID- 6)  50,  50,  bO 
50  WRIIE(o,1G0)X,Y(1),Y(2),Y(3),¥(4)  ,Y(5) 

100  FORMAT ( 1X,6C20.6) 

XP=XF+Q . 05D00 
60  DO  20  1=1,5 
20  Y  (I)  =  Y  (I)  ♦H*DERY  (I) 

X  =  X  >  h 

IF(X. GT.8.00D00) GO  TO  30 
GO  TO  10 
30  STOP 
END 


nr.  n  n  nc.  nnnnr.  nn  n  n 


1  IV  G  LEVEL  21 


MAIM 


DATE  =  80  1 89 


14/30/26 


C 

SUBROUTINE  CAC(X,Y,Z) 
IMPLICIT  REAI.  +  8  (A-H,P-Z) 
DIMENSION  Y  (5)  ,Z  (5,6) 
COMMON  GC, AR,CK, AM, S 


IF(X.GT. I.OODOO) GO  TO  11 
E  =  1  »  OODOO *0 . 50D0  C’*X 
GO  TO  12 

11  B  =  1 .50D00*(1. 00D00-0. 1 0DOO  *  X ) 


12  A— 0 . 20D00*X ♦  1 . 00  D00 
D  A- 0.20 DO 0 


U— 0. 12D00*X*0.  30D-3 
DU=  0.12D00 


If  (J.LT.  0.30000)  GO  TO  31 

U=0.30D00 

DU=  0,00100 

31  CSP=DSQRT  (Y  (5) ) 

S IH= (Y ( 4) **10) /DSP 
SIC= (Y (3) **10) /DSP 
5TH=DSQBT  (Y  (4) ) 

STC=DSQRT (Y (3) ) 

EETAW  =  B*STH/Y  (5) 

EET  AC=B*STC/Y (5) 

ARK=AR*CK 

AA1=SIH/(  1. OODOO ♦ BETA W*bET A  W) 

AA2=SIC/ (1. OODOO ♦BETAC*BET AC) 

AKC=AA1* (CK*Y (1) -Y (2) ♦ ARK* DETAN *Y ( 1 ) )/Y (1 ) /A A2+ 1.0  0DO0-A3K*uETAC 
AA3=ARK*ARK+ (AKC- 1. OODOO) *  (AKC-  1.00D00) 

AA=Y  (1)  *D*U 
GG*GC/(GC-1, OODOO) 

CC=GC*AM*AM 
Z  (1  ,  1)  =0.  OODOO 
Z  (1,2)  *0.  OODOO 
Z  (1,J)=GG*Y  (5)/Y  (3) 

Z  (1 ,4)*0.  OODOO 
Z  (1 ,5)  =-1. OODOO 
Z  (1,6)=S*AA2*AA3*AA 


AA4=ARK*BETAC>AKC-1,0UD00 
Z  (2,1)  *CC*Y  (5)  *Y  (1)/Y(3) 

Z  (2,2)  =0.  OODOO 
Z  (2,3)  =0.  OODOO 
Z  (2,4)  =0.  OODOO 
Z  (2,5)*  1.  OODOO 


U  IV  G  LEVEL 


DATE  =  60  1  69 


14/30/26 


2  (2,6)=S*AA2*AA4*AA 


AA5=CC*  (Y  (1)  -Y  (2)  )  *  (Y  ( 1 )  -Y  (2) )  *0.  5OD00-GG*  (Y  (4)  -Y(3)) 
AA6=Aa:>*ARK*Y(1)  *Y  (  1)  ♦  (CK*Y  ( 1)  -  Y  (2)  |  ♦  (CK*Y  (1 )  -Y  (2)  ) 
66=6*3 

Z(3,1)  =  0.00D00 
Z  (3,2) =0. OODOO 
Z  (3 ,3)  =0 . 00 DOO 
Z  (3,4)  =U*GG*Y  (5)  *Y  (2)  /Y  (4) 

Z  (3,5)  =-0*Y  (2) 

Z  (3,6)  =Y  (5)  *AA5*DU*Y  (2) /Y  (4)  +U*6*AA1*AA6*B6 

AA7=CK*Y  (1)  *  (AP* BETAS*  1-00D00)  ”Y(2) 

FP=Y(5)*Y(2)/Y(4) 

Z  (4 , 1)  =  0.  OODOO 
Z(4,2)=0*CC*PP 
Z  (4,3)  =0. OODOO 
Z  (4 ,4)  =0.  OODOO 
Z  (4  5)  =0 

Z (4*6) =CC*PF*DU*  ( Y  ( 1 ) “Y  (2) ) ♦  U* S* A A  1  * A A7*B6 


AA8= (Y ( 1) *Y  ( 4 )  “ Y  (2) *Y  (3) ) 

AA9=Y(2)  *Y(3)  *U*AA6 

Z(5,1)  =  (AA9-0*Y(1)*Y(4)  )/(Y  (1)*AA9) 

Z  (5,2)  =  (AA9-Y  (2)  *Y  ( 3)  ♦  ( 1 . 00  DOO-U)  )  /  (Y  (2)  *  A  A9 ) 
Z  (5,3)  =-Y  (2)  *(1.0OD00-U)/AA9 
Z(5,4)=-Y  (1)  *U/AA9 
Z (5,5) = 1 . 00  DOO/Y  (5) 

Z  (5,6)  =AA8*DU/AA9-DA/A 

FETUBN 

END 


i  IV 


LEVEL  21 


MAIN 


DATE  =  80189 


14/30/26 


?6 


c 

c 

SUBROUTINE  FCT  (X  ,Y , DEFY) 
IMPLICIT  KEAL*8  {A-H,P-Z) 
DIMENSION  Y  (5)  ,0  ERY  (5)  ,  A (5 , 6) 
CCHHON  GC,AR,CK, ftd,S 
C 
C 

CALL  CAC  (X,  Y, A) 

N=5 
M=N  ♦  1 
C 
C 

DO  30  K=1,N 
KFl=K+1 

IF (K.EQ. N) GO  TO  11 
J  J=  K 

EIG=DABS (A  (K, K) ) 

DC  7  I  =  KP  1 , N 
AB=  DABS  (A  (I ,  K)  ) 

IF(BIG-AB)  3,7/7 
3  E JG=AB 
JJ=I 

7  CCNIINOE 
IF(dJ-K) 8,1 1.8 

8  DO  9  J=  K,  d 
IEMP— A  (JJ  ,J) 

A (JJ.J) =A  (K,J) 

9  A(K,J)=TEttP 

11  DO  20  J=KP1 «M 

20  A  (K  ,  J)  =A(K,J)/A  (K,K) 

DO  30  1=1, N 
IF(I.EU-K)GO  TO  30 
CG  25  J=KP1 ,d 

25  A(I,J)SA(I,J)~A(I,K)*A(K,J) 

30  CONTINUE 
DC  40  1=1,  N 
40  DER  Y  ( I)  =  A  (I  ,  d) 

RETURN 

END 
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LIST  OF  SYMBOLS 
(except  Section  1.2) 


2 

Cross-sectional  area  of  channel*  m 

2 

Cross-sectional  area  of  core  region,  m 

2 

Cross-sectional  area  of  wall  layer,  m 
Magnetic  field  strength,  T(esla) 

Width  of  generator  duct,  m 
Specific  heat  capacity,  J4kg*K) 

Electric  field  strength,  N/C 
CpT  +  u2/2,  J/kg 
Specific  enthalpy,  J/kg 

2 

Electric  current  density,  C/(m  *s) 

Loading  factor 

VucB 

Ey*/uwB 
Mach  number 

Total  mass  flow  rate  through  generator,  kg/s 

Mass  flow  rate  In  core  region,  kg/s 

Mass  flow  rate  In  wall  layer,  kg/s 

Unit  normal  vector 
2 

Pressure,  N/m 

Gas  constant,  J/(kg-K) 

W10  Interaction  parameter 
Specific  entropy,  J/kg*K 
Temperature,  K 


?<? 


U  Specific  Internal  energy,  J/kg 

u  Flow  velocity,  m/s 

V  Specific  volume,  m3/kg 

x  Cartesian  coordinate  along  the  duct  axis,  In  the  flow  direction, 

measured  from  the  entrance,  m 

o  Ratio  of  Hall  field  to  Faraday  field  (E  /E  ) 

x  y 

e  Hall-current  parameter 

y  Ratio  of  specific  heat  capacities  (cp/cv) 

$  Angle  of  Inclination  of  segments  (=  tan”1aw.  Figs.  2.1  to  2.2),  rad. 

u  Ratio  of  mass  flow  rate  In  wall  layer  to  the  total  mass  flow  rate 

through  the  generator 

v  Ratio  of  cross-sectional  area  of  wall  layer  to  that  of  entire  duct 

3 

p  Density,  kg/m 

-1  3 

a  Electric  conductivity,  n  /m 

Superscripts 

(  )'  Dimensionless  quantities  (see  Section  3.2  and  beginning  of  Section 
3.1) 

Vectorial  quantities 

A 

(  )  Unit  vectors 

Subscripts 

c  Of  core  region 

1  At  entrance 

w  Of  wall  layer 

x  Component  in  the  x-dl recti on 

y  Component  In  the  y-dlrection  (see  Fig.  2.1) 


I 


